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FOREWORD 


The  work  described  in  this  report  wns  authorized  under 
Task  1 B522301AQ8101,  Di  .-semination  Investigations  of  Liquid 
and  Solid  Agents  (U).  The  work  was  started  in  July  1966  and 
completed  in  January  1967. 

Reproduction  of  this  document  in  whole  or  in  part  is  pro¬ 
hibited  except  with  permission  of  the  CO,  Edgewood  Arsenal, 
ATTN:  SMUEA-RPR,  Edgewood  Arsenal,  Maryland  21010;  how¬ 
ever,  Defense  Documentation  Center  is  authorized  to  reproduce 
the  document  for  US  Government  purposes. 

The  information  in  this  document  has  not  been  cleared 
for  release  to  the  general  public. 


DIGEST 


The  conditions  required  for  the  spontaneous  ignition  of 
a  volatile  fuel  drop  suddenly  exposed  to  a  hot  oxidizing  en¬ 
vironment  have  been  investigated  theoretically  using  a  simpli¬ 
fied  model  of  the  physical  situation.  Analysis  of  the  model, 
a  reactive  spherical  shell  with  unequal  boundary  temperatures, 
defined  the  conditions  consistent  with  thermal  stability  of 
the  shell.  Two  types  of  reaction  rate  were  used,  one  with  a 
simple  Arrhenius  form,  and  a  second  explicitly  involving 
radial  distance  as  well  as  a  temperature  dependence  in 
Arrhenius  form. 
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I  INTRODUCTION 


Many  physical  and  chemical  processes  are  involved  in  the  ignition 
of  a  volatile  fuel  droplet  suddenly  exposed  to  a  hot  oxidizing  environ¬ 
ment.  A  complete  mathematical  description  of  these  processes  and  their 
interactions  can  be  formulated  but  is  so  complex  that  even  computer 
solutions  are  impractical.  Consequently,  simplification  is  desirable. 

A  simplifying  approach  to  the  problem  is  suggested  by  an  experimental 
study  of  droplet  ignition. 1  That  study  revealed  that  igniti-n  only 
occurs  if  the  initial  droplet  diameter,  d i ,  exceeds  a  critical  diameter, 
d  ,  dependent  on  the  nature  of  the  fuel  and  on  the  oxidizer  temperature 
In  the  case  of  hexadecane  droplets,  d  decreases  with  increasing  oxidiz 
temperature  but  is  insensitive  to  oxidizer  concentration  for  concentra¬ 
tions  greater  than  the  minimum  required  for  flammability.  These  obser¬ 
vations,  in  particular  the  existence  of  a  critical  initial  diameter, 
suggest  a  similarity  between  drop  ignition  and  smmpler  steady-state 
problems  involving  attainment  of  a  critical  state  in  which  heat  released 
by  chemical  reaction  cannot  be  dissipated  to  the  environment. 

Analyses  of  the  thermal  stability  of  reactive  materials  in  the  form 
of  slabs,2  cylinders,  and  spheres3  involved  a  criticality  parameter  S 
related  to  the  problem  parameters  by  Eq.  1: 

o  =  -  l 2  q  0  ( 1  ) 

MTl 

whe  re 

E  -  an  overall  activation  energy 
Ft  =  universal  gas  constant 
Tq  =  boundary  temperature 

q0  -  rate  of  heat  release  per  unit  volume  at  70 
(in  Arrhenius  form) 

l  -  a  characteristic  dimension  of  the  sample 
K  -  sample  thermal  conductivity. 


If  b  exceeds  a  certain  critical  value,  §cr,  dependent  on  sample  geometry, 
the  material  will  ignite  or  explode  in  a  time  that  is  related  to  q0 ■  for 
a  slab  with  boundaries  at  different  temperatures,  q 0  refers  to  the  reactjo 
rate  at  the  hot  boundary  and  8 f r  depends  on  the  temperature  difference  be¬ 
tween  the  two  boundaries.4 

By  analogy  to  the  cited  cases  the  critical  state  for  droplet  ignition 
will  involve  at  the  least:  (1)  the  surface  temperature  of  the  droplet, 

(2)  the  ambient  temperature,  (3)  the  size  of  the  droplet,  and  (4)  the 
kinetics  of  the  reaction  between  evaporated  fuel  and  gaseous  oxidizer. 
Because  the  rate  of  reaction  will  be  affected  by  the  concentration  gra¬ 
dients  associated  with  the  movement  of  evaporated  fuel  away  from  the 
drop,  the  evaporation  rate  will  also  be  involved  if  only  indirectly.  A 
characteristic  dimension  analogous  to  the  term  l  in  Eq .  1  is  not  well 
defined  for  a  process  involving  diffusion  from  a  finite  source  into  an 
infinite  medium.  An  artificial  outer  boundary  must  be  imposed  in  order 
to  prescribe  definite  boundary  conditions. 

This  theoretical  study  is  concerned  with  two  simplified  models  of 
droplet  ignition.  Both  models  involve  a  steady-state  approach  (in  which 
time  does  not  appear  as  a  variable)  to  the  determination  of  the  thermal 
stability  of  a  reactive  spherical  shell  with  an  outer  boundary  temperature 
greater  than  that  of  the  inner  boundary.  The  first  model,  the  simpler 
case,  involves  a  reaction  rate  in  Arrhenius  form  that  varies  explicitly 
with  temperature  but  not  with  position  except  as  temperature  varies  with 
position.  The  second  model,  a  more  realistic  and  more  complicated  one, 
includes  a  reaction  rate  that  varies  explicitly  with  both  position  and 
temperature.  The  form  of  the  dependence  on  position  was  chosen  to  simu¬ 
late  the  variation  of  fuel  concentration  with  distance  from  an  evaporating 
drop.  In  both  cases,  the  intention  was  to  determine  8  r  as  a  function  of 
the  various  parameters  which  enter  the  problem. 
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II  DISCUSSION 


A.  Model  I:  Reaction  Rate,  a  Function  of  Temperature 

For  a  steady-state  temperature  distribution  to  exist  in  the  situation 
represented  by  Fig.  1,  the  temperature  within  the  shell  must  satisfy  the 

differential  equation 

W2T  =  -g  (2) 

where  \  is  the  thermal  conductivity  of  the  shell  material  and  a  is  the 
rate  of  heat  release  per  unit  volume  within  the  shell.  The  prescribed 
boundary  conditions  are  T  -  7  at  r  -  a,  and  7  =  T,  at  r  =  b.  If  the 
rate  of  heat  release  q  is  assumed  t<o  have  the  Arrhenius  form  q  -  Ac~“  Kr, 
a  steady -state  temperature  distribution  is  not  possible  for  all  values  of 
the  parameters  u ,  b,  T a ,  7 b ,  A,  E,  and  K.  The  task  is  to  determine  the 
range  of  values  that  do  permit  a  steady-state  temperature  distribution. 


TA  -  ©00 


FIG.  1  REACTIVE  SPHERICAL  SHELL 
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Solutions  of  the  System  of  Eqs.  (6)-(8)  are  only  possible  for  values  of 
o  less  than  a  critical  value  8  that  depends  on  the  boundary  conditions. 
Because  analytic  solution  of  Eq.  6  is  not  possible,  specific  values  of 
0  must  be  obtained  bv  numerical  methods.  These  values  of  ?  define 

c  r  '  c  r 

the  possible  ranges  of  variation  of  !  ,e  parameters  a,  b,  T  ,  Tk,  A,  E, 
and  K . 
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To  determine  the  maximum  value  of  5,  8fr  =  8  c  .  (  z  0 ,  t~0 )  which  permits 
the  solution  of  differential  Eq.  (6),  subject  to  the  boundary  conditions 
(7)  and  (8),  the  following  method  was  adopted: 

Starting  with  a  definite  value  of  5,  Eq.  (6)  was  integrated  subject 
to  the  boundary  conditions 

6  =  0  and  — —  =  m  at  y  =  0  (10) 

dy 

where  m  is  a  parameter.  The  solution  6  =  9(y,S,m)  is  a  function  of 

the  variable  y  and  the  parameters  S  and  m.  The  value  of  6  at  y  -  y0  is 
0(y  8,  m)  .  As  m  varies,  this  quantity  has  a  maximum  value  8  -  0(yo,5), 

i .  e  .  , 


~Q  =  f(v  g)  =  max  9(y0,S,m) 

■yCI  '  -CD<*  <CO 


(11) 


If  (9  is  less  than  6Q,  it  is  evident  that  Eq.  (6),  subject  to  the  boundary 
conditions  (7)  and  (8),  has  no  solution  for  the  particular  value  of  § 
chosen,  i . e . ,  8  <  r (^0 1 '  If  8  is  greater  than  9 ^ ,  there  are  two 

solutions  of  the  system  of  Eqs.  ( 6 ) - ( 8 )  and  S  <  §cr-  •  Finally,  if  9  is 
equal  to  9^,  there  is  only  one  solution  to  the  problem  and  the  value  of 
8  is  the  critical  value  S g r  corresponding  to  given  values  of  9 Q  and  y0 
(or  r  0 ) - 

As  noted,  the  numerical  procedure  provides  two  solutions  to  Eq.  6 

if  S  <  S  Other  studies6'?  have  revealed  multiple  solutions  in  cases 

c  r 

involving  spherical  symmetry.  In  particular,  a  large  number  of  solutions 
(not  necessarily  stable)  are  possible  for  some  values  of  8.  We  are  con¬ 
cerned  in  this  study  with  values  of  8  in  the  vicinity  of  Scr  for  which 
there  are  probably  at  most  two  solutions.  Of  these  two  solutions  the  one 
corresponding  to  the  smaller  value  of  m  is  probably  the  significant  one, 
for  in  a  physical  problem  the  stable  temperature  distribution  would  be 
approached  from  below  and  the  eventual  steady-state  distribution  would 
be  that  corresponding  to  the  lowest  value  of  n. 

A  computer  program  based  on  the  method  described  was  carried  ojt  to 
determine  5fr  as  a  function  of  zQ  and  6>0 .  The  results  are  shown  in  Table  I. 
The  range  of  values  chosen  for  9 g  (0  to  -14)  corresponds  to  the  range  of 
values  for  T  TL.  and  E  that  might  be  encountered  in  an  experimental  study 
of  drop  ignition.  Graphical  interpolation  may  be  used  to  obtain  other 
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Table  I 


CRITICAL  VALUES  OF  8:  MODEL  I 


<£> 

O 

II 

o 

0O  =  -2.0 

e0  =  -4.0 

*0  =  -6-°j 

&0  =  “8.0 

o 

II 

l 

o 

o 

60  =  -14. ( 

=  0.  05 

|  3.74 

4.21 

4.69 

5.16 

5.64 

6.13 

7.  13 

=  0.125 

4.47 

5.72 

7.04 

8.43 

9.90 

11.45 

14.79 

=  0.25 

6.16 

9.21 

12.75 

16.75 

21.21 

26.12 

37.30 

=  0.5 

14.  00 

26.01 

42.30 

62.75 

87.31 

115.92 

185. 25 

=  0.75 

56.18 

121.78 

222.66 

360.04 

534.00 

744.20 

1273.2 

=  0.90 

351.34 

820.45 

1591.2 

2685.5 

4107. 1 

5855.3 

10328.5 

values  of  8  for  8  from  0  to  -14  and  for  zn  from  0.05  to  0.90.  Ex- 
trapolation  of  the  tabulated  results  to  values  of  zQ  in  the  range 
0  £  zQ  <  0.05  is  e Iso  readily  accomplished.  As  suggested  by  the  data  in 
Tables  I  and  II,  for  all  (?0 ,  as  zQ  -*  0,  Scr  appears  to  approach  3.32,  the 
value  for  a  solid  sphere.3  The  convergence  was  not  investigated  analyt¬ 
ically.  Extrapolation  of  the  data  to  values  of  zQ  >  0.9  is  not  feasible. 
However,  for  1  >  z0  >  0.9,  the  reactive  shell  of  Fig.  1  may  be  approxi¬ 
mated  by  a  thin  slab  and  the  solution  to  the  thin  slab  problem4  used  to 
estimate  values  of  8  for  the  shell.  To  utilize  the  slab  solution  it 

c  r 

must  be  noted  that  8  for  a  slab  is  defined  in  terms  of  the  characteristic 
dimension  (6  -  a) / 2 . 


Table  II 

CRITICAL  VALUES  OF  5:  MODEL  I 


*0  =  0  90 


e0  = -2 


zn  =  0.001  3.330  |  3.335  3.340 


=  0.01  3.403 


*0  =  0.02  3.485 


n 


B.  Model  I:  Drop  Ignition 


In  order  to  apply  Model  I  criticality  data  to  drop  ignition  it  is 
necessary  to  consider  the  variation  of  a,  the  drop  radius,  and  of  b,  the 
radius  of  the  effective  outer  boundary,  with  time.  The  surface  tempera¬ 
ture  of  a  drop  suddenly  exposed  to  a  hot  oxidizing  atmosphere  will  quickly 
approach  a  value  Ta  corresponding  to  a  balance  between  heat  transfer  in¬ 
ward  and  diffusion  of  evaporated  fuel  outward.  Thereafter,  the  surface 
temperature  Ta  will  remain  nearly  constant  as  evaporation  proceeds  and  the 
radius  of  the  drop  decreases.  Initially,  the  radius  b  of  the  effective 
outer  boundary  equals  a  and  then  increases  as  a  decreases  and  evaporated 
fuel  diffuses  into  the  surrounding  atmosphere.  The  outer  radius  b,  although 
poorly  defined,  must  eventually  approach  a  maximum  value  corresponding  to 
completr  evaporation.  As  shown  schematically  in  Fig.  2,  during  this  time, 

5  will  approach  a  maximum  value  as  a  approaches  zero.  If  this  maximum 


TA  -  4900-  #01 


FIG.  2  DROP  TRAJECTORY  IN  S  -  Z0  PLANE:  MODEL  I 
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value  exceeds  3.32,  the  apparent  limit  of  3  as  z 0  0  for  all  6^,  the 

drop  can  ignite.  For  fixed  Tfc  ,  the  above  qualitative  description  implies 
the  existence  of  a  d  r  such  that  a  drop  with  an  initial  diameter  dt  <  dcr 
will  not  ignite. 

As  noted,  all  3  r  as  a  function  of  6 0  appear  to  approach  the  value 
of  80  =  3.32  as  z0  -  0.  If  the  final  value  of  6  corresponding  to  complete 
evaporation  of  a  drop  is  proportional  to  the  initial  drop  diameter  dx,  then 
the  definition  of  3  implies  that  dcr  will  decrease  as  Tb  increases  with  an 
apparent  activation  energy  that  can  be  deduced  from  the  definition  of  o 
(Eq.  9).  That  definition  results  in  the  proportionality: 


KRTl 


E /  RT  . 


(12) 


The  pre-exponential  term  A  is  not  strictly  independent  of  temperature  be¬ 
cause  it  includes  an  average  reactant  concentration  that  varies  '-.lth  T a. 
If  Th  »  Ta,  the  dependence  of  A  on  Tb  is  slight  and  can  be  ignored  with 
the  result  that 


d  In  d c  r 

d(l/Th) 


(13) 


where  7\  is  a  mean  oxidizer  temperature  in  the  temperature  range  of  in¬ 
terest  and  Q  is  an  apparent  activation  energy.  The  measured  value  of  Q 
for  ignition  of  hexadecane  drops  in  hot  air  is  7  kcal/mole,7  corresponding 
to  a  value  for  E  of  15  kcal/mole. 


C,  Model  II:  Reaction  Rate,  a  Function  of  Position  and  Temperature 

Analysis  of  Model  I  revealed  the  dependence  of  5  p  on  the  dimension¬ 
less  temperature  difference  between  the  hot  and  cold  boundaries  of  the 
reactive  spherical  shell.  The  study  also  revealed  that  the  Spp  of  impor¬ 
tance  in  determining  the  drop  critical  initial  size  was  that  corresponding 
to  a  solid  sphere,  Scr  =  3.32.  The  model,  however,  is  unreal  because  it 
does  not  include  the  effect  of  concentration  gradients  associated  with  the 
diffusion  of  fuel  away  from  an  evaporating  drop.  If  the  reaction  rate  q 
varies  explicitly  with  fuel  concentration  ,  with  oxidizer  concentration 
C0  ,  as  well  as  with  temperature  T,  the  concentration  gradients  will  affect 
the  conditions  required  for  criticality. 
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The  effect  of  variable  reactant  concentration  can  be  simulated  by 
including  in  the  model  a  rate  of  heat  release  given  by 


q  =  ACfC0xe' 


In  the  physical  situations  of  interest  the  dimensionless  boundary  tempera¬ 
ture  <<  0.  By  analogy  to  the  slab  problem,4  for  dQ  <<  0,  the  maximum 
reaction  rate  occurs  in  the  vicinity  of  the  hot  boundary.  In  that  region 
the  concentration  of  oxidizer  is  approximately  equal  to  the  value  in  the 
ambient  atmosphere  and  can  be  absorbed  into  tne  p re - expon en t i a  1  term  A, 


giving 


where 


AlCfe~E/RT 


The  reaction  rate  in  Model  II  is  assumed  to  be  of  the  form  given  in 
Eq.  (15).  The  boundary  conditions  are  T  -  Ta  at  r  =  a  and  T  =  Tb  at  r  =  6. 
The  element  that  differentiates  Model  II  from  Model  I  is  the  dependence  of 
the  reaction  rate  on  position  through  the  fuel  concentration  Cj.  Because 
the  variation  of  Cj  with  position  is  intended  to  simulate  the  results  of 
diffusion,  the  concentration  Cj  is  required  to  satisfy  the  equation 
V2C^.  =  0  with  the  boundary  conditions  =  C0  at  r  -  a  and  =  0  at 
r  =  6 .  Tha  t  is, 


.ij 


In  terms  of  the  variable  y,  Eq .  16  becomes 


yd  -  y0) 

y  ‘  (1  -  y) 


where  yQ  =  1  -  (a/6)  =  1  -  zQ.  Now  let 


w; 
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Then  Eq.  2  becomes 


d2e  2  dd 
dy2  1  '  y  dy 

with  boundary  conditions 

0=0  at  y  =  0 

0  *  0O  <  0  at  y  =  y Q  (20) 

As  before,  the  problem  is  to  determine  as  a  function  of  and  z 0  the 
maximum  value  of  S,  8 c  ,  that  permits  solution  of  Eq.  19.  The  range  of 
variation  of  6g  and  z0  is  the  same  as  in  the  analysis  of  Model  I. 

Determination  of  specific  values  of  8fr  was  done  in  the  same  manner 
as  in  the  analysis  of  Model  I  with  the  results  shown  in  Table  III..  Com¬ 
parison  of  the  data  in  Tables  I  and  III  reveals  that  the  inclusion  of  a 
concentration  gradient  has  resulted  in:  (1)  a  marked  increase  in  Scr  for 
fixed  0Q  and  zQ,  (2)  a  minimum  in  8  r  as  a  function  of  z  0  for  fixed  0g  , 
and  (3)  a  rapid  increase  in  8  ^  r  as  zQ  decreases  towards  zero  from  the  value 
corresponding  to  the  minimum  in  S  .  The  cited  minimum  value  is  quite 
evident  for  small  values  of  0g.  However,  it  occurs  at  a  value  of  zQ  that 
decreases  with  decreasing  8g,  and  therefore  it  is  not  apparent  in  Table  III 
for  large  | 9 g  |  .  The  rapid  increase  in  §cr  with  decreasing  z0  for  zQ  less 
than  the  value  corresponding,  to  the  minimum  in  8  was  not  investigated 
in  detail. 


Table  III 

CRITICAL  VALUES  OF  8:  Model  II 


o 

II 

o 

0„=  -2 

e0=  -4 

£„=  "6 

OO 

1 

It 

o 

V  -I0 

1 

II 

O 

zQ  =  0.05 

50.48 

66.29 

83.01 

100.80 

119.76 

139.96 

184.32 

z0  =  0.125 

29.09 

46.06 

66.33 

91 . 52 

120.33 

153.46 

233.65 

,  -  0.25 

‘0 

24.73 

47.27 

79.77 

123.40 

179.35 

248.87 

433.60 

z.q  =  0.50 

37.35 

L. 

89.26 

181.88 

327.61 

538.1  9 

825.10 

1673.72 

z q  =  0.75 

122.62 

337.63 

781 .56 

1562.96 

2788.89 

4563.9 

10174.7 

z0  =  0.90 

707.39 

2078.43 

5114.69 

10764.0 

19992.4 

33746.3 

78559.5 

-8e 


yd  -  V 
y0(1  -  y) 


(19) 
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D.  Model  II:  Drop  Ignition 

The  existence  of  a  minimum  §cr  for  fixed  80  affects  the  initial  drop¬ 
let  size  required  for  ignition.  A  drop  trajectory  in  the  S-zQ  plane  must 
intercept  the  8  -z^  curve  for  ignition  to  occur.  As  shown  in  Fig.  3, 

the  critical  initial  droplet  size  is  that  size  corresponding  to  the  tra¬ 
jectory  tangent  to  the  8^r-z0  curve.  Because,  in  practical  cases,  8Q  <<  0, 
the  point  of  tangency  will  occur  at  values  of  z  only  slightly  greater  than 
zero.  Unlike  the  situation  represented  by  Model  I,  the  §cr  of  interest 
is  greater  than  3.32  and  increases  with  decreasing  8Q. 

Because  the'S  that  determines  drop  critical  size  varies  with  9^ 
and  occurs  at  values  of  z Q  >  0,  deduction  of  the  temperature  variation 
of  d  is  more  difficult  than  it  was  in  the  case  of  Model  I.  However, 

e  r 

in  the  cases  of  interest,  dQ  «  0  and  the  critical  droplet  trajectory  is 
tangent  to  the  criticality  curve  at  a  value  of  zQ  very  close  to  zero. 
Consequently,  it  may  be  assumed  that  the  value  of  b  at  the  point  of  tan¬ 
gency  is  proportional  to  d^  r  .  The  definition  of  8,  Eq.  9,  then  leads  to 


(21) 


where  8 T ,  the  value  of  8 c r  at  the  point  of  tangency,  will  be  very  close 
to  the  minimum  value  of  8  for  a  given  6n.  If  T  ,  Cn  ,  and  Cn  are  held 

cr  °  0  a  Ox  U 

constant,  then  the  apparent  activation  energy  Q  is  given  by 


a  relation  differing  from  Eq .  13  in  the  third  term  on  the  right.  That 
term  can  be  re-expressed  in  the  form 


IS 


*0 

TA-  <*00  -  00? 


FIG.  3  DROP  TRAJECTORY  IN  5  -  Z0  PLANE:  MODEL  II,  0O  =  ~6 

The  derivative  of  60,  with  Ta  constant,  gives 


and  may  be  either  positive  or  negative  depending  on  the  specific  values 
of  Ta  and  7’ft.  From  the  data  in  Table  III  for  values  of  that  show  a 
minimum  in  the  §rr-z0  curve  a  value  for  the  term  d  In  5 T/ddQ  can  be  esti¬ 
mated  as  about  5  *  10  2.  Consequently,  the  contribution  of  the  third 
right-hand  term  in  Eq.  22  to  the  overall  activation  energy  is  approximately 


In  many  cases  of  interest,  (Tt  -  Ta ) IT h  =  1  and 

01  =  -5  x  10'2  —  (26) 

2 

In  such  cases  Q'  makes  only  a  small  contribution  to  Q  and  can  be  ignored. 
In  other  cases,  the  absolute  value  of  |2(T0/Tfc)  ~  1|  may  be  substantially 
less  than  one  and  the  contribution  to  Q  will  be  even  less  important.  The 
result  is  that  both  Models  I  and  II  lead  to  a  value 

Q  =  £  -  RTb  .  (27) 


1 1 1  SUMMARY 


The  results  of  the  analyses  of  both  Models  1  and  II,  when  applied  to 
drop  ignition,  predict  the  existence  of  a  critical  initial  diameter,  d  , 
such  that  a  drop  with  an  initial  diameter  di  >  dcr  will  ignite  prior  to 
complete  evaporation,  whereas  a  drop  with  dL  <  dcr  will  merely  evaporate. 
For  equivalent  conditions,  d^r  for  Model  I  will  be  substantially  less  tiian 
d  for  Model  1 1 .  Both  models  predict  that  d  will  decrease  with  increas- 
ing  7\  with  an  overall  activation  energy  Q  approximately  equal  to 


Model  I,  while  formally  interesting,  does  not  include  the  concentration 
gradients  that  are  an  inherent  part  of  the  phvsical  problem.  Consequently, 
the  usefulness  of  the  Model  I  analysis  lies  primarily  in  determination  of 
the  dependence  of  on  9 Q  and  ZQ .  Model  II  includes  a  variation  of  fuel 

concentration  with  distance  from  the  source  and  therefore  more  closely 
resembles  the  physical  situation.  The  analysis,  however,  is  not  general 
because  of  the  form  assumed  for  the  reaction  rate.  Inclusion  of  a  fuel 
concentration  variation  has  two  effects.  The  criticality  parameter  S,  as 
defined  in  F.q.  18,  now  explicitly  depends  on  CQx  and  CQ  .  Because  C0  x  and 
CQ  appear  in  the  p re  -  exponent ia  1  factor,  the  effect  is  minor  and  similar 
for  both  fuel  and  oxidizer.  A  variation  in  either  C„  or  C„  can  be  nul- 
lified  by  a  much  smaller  variation  in  Tb.  The  second  effect  is  more  im¬ 
portant.  The  pronounced  increase  in  5  f  resulting  from  the  variation  of 
Cj  with  distance  from  the  drop  is  an  indirect  consequence  of  using  a  reac¬ 
tion  rate  proportional  to  .  A  reaction  rate  proportional  to  a  power  of 
>  1  would  presumably  enhance  the  effect.  Consequently,  the  reaction 
order  with  respect  to  fuel  is  probably  an  important  parameter  in  deter¬ 
mining  critical  drop  size.  The  importance  of  the  reaction  order  with 
respect  to  oxidizer  is  much  less  because  the  maximum  reaction  rate  occurs 
near  the  hot  boundary  where  the  concentration  of  oxidizer  is  nearly  the 
same  as  that  in  the  ambient  atmosphere  and  is  not  sensitive  to  diffusion. 
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1  3  A05TRAC  T 


The  conditions  required  for  the  spontaneous  ignition  of  a  volatile  fuel 
drop  suddenly  exposed  to  a  hot  oxidizing  environment  have  been  Investigated 
theoretically  using  a  simplified  model  of  the  physical  situation.  Analysis  of 
the  model,  a  reactive  spherical  shell  with  unequal  boundary  temperatures,  de¬ 
fined  the  conditions  consistent  with  thermal  stability  of  the  shell.  Two 
types  of  reaction  rate  were  used,  one  with  a  simple  Arrhenius  form,  and  a 
second  explicitly  involving  radial  distance  as  veil  as  a  temperature  dependence 
in  Arrhenius  form. 
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